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Abstract 
The B — «£v weak decay form factors are calculated via light-cone sum rules within the frame- 
work of the heavy quark effective theory. We calculate the leading and the relevant sub-leading 
universal form factors. Our results are matched to the known soft pion limit. We also address 
the large pion energy limit of our sum rule results. Our results are compared with that of other 
approaches. 
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I. INTRODUCTION 


It is well-known that the weak B — 7 transition is important for the extraction of 
the Cabbibo-Kobayashi-Maskawa matrix element V,,, from semileptonic decays, and for the 
measurement of CP violation from non-leptonic decays. In order to achieve these goals 
advances have to be made both on the experimental and theoretical sides. On the theoretical 
side an important task is to reduce the uncertainties in the calculations of the relevant 
hadronic matrix element, represented by 


2 m2 2. a2 
(x(p)| a b|B(P)) = f(g?) Pia sg + fol?) a (1) 


where q = P — p is the momentum transfer to the leptons. Up to now, the form factors 
have mostly been calculated via light-cone (LC) sum rules [I] in full QCD. Instead, here we 
will work in the context of the heavy quark effective theory (HQET) [B]. The reasons are as 
follows: 

(i) Unlike for heavy-to-heavy transitions an application of HQET methods does not sig- 
nificantly simplify the analysis for heavy-to-light transitions. Nevertheless, the systematic 
nature of HQET allows one to identify and estimate uncertainties in the heavy-to-light 
transitions more easily. 

(ii) For an analysis involving the B meson HQET is the correct approximation method 
of QCD. In certain parts of the phase space of the semileptonic decay, the pion is not very 
energetic so that HQET may still be valid. 

(iii) The B — mm amplitude has been shown to be factorizable at the leading order of 
the 1/m» expansion [B]. Knowledge of the form factors based on HQET is needed so as to 
consistently apply the factorization approach. 

(iv) Furthermore, it is still controversial whether the time-like transition B — 7 at large 
recoil is governed by perturbative or non-perturbative QCD |p, HJ. If non-perturbative effects 
dominate, then the use of HQET for the B — 7 transitions can be fully justified. 

(v) In addition, the results of HQET can also be applied to D — «(K) transitions. In 
HQET, a model-independent analysis of B(D) — 7 transitions to order 1/mg has been 
done in Ref. P, B] and for B(D) — p in Ref. [7]. 

HQET simplifies the analysis by introducing a set of universal functions. However, in 
order to obtain information on the universal functions themselves, some nonperturbative 
techniques, such as light-cone sum rules or lattice simulations, must be used in addition. 
Here we adopt the LC sum rule method fÑ, B] which is suitable for the calculation of form 
factors when light energetic hadrons are involved. To our knowledge, the LC-HQET sum 
rules were first applied in Ref. [10]. In Ref. [I1], they were used at the leading order of 
HQET for B — - transitions. However, the results in this paper differ from those in [11]. 

Note that we distinguish between LC-HQET sum rules and LC-QCD sum rules with the 
1/mg expansion [12]. The main reason is that the ways to include radiative corrections are 
different. Another reason is that for H — 7 transitions, there is a subtle difference between 
the two types of sum rules, as will be discussed in the paper. 

In this paper we apply LC-HQET sum rules to the H — 7 transition to order 1/mg. In 
the next section we give a brief review on the application of HQET to the decay B — «4v. In 
section [I], the leading and next-to-leading order universal functions are calculated by using 
LC-HQET sum rules. In Section [V| we present our numerical results. Section [V] contains a 
summary and discussion, and a comparison with other approaches. 


Il. THE HEAVY QUARK EXPANSION 


In HQET, the velocity of heavy quark Q, v, is a well defined quantity. The heavy quark 
field can be represented by the velocity-dependent field, 


hy(x) = exp(img v-x) P, Q(x), (2) 


1 
where P, = LU projects onto the upper component of the heavy quark field Q(x). To 
the order 1/mg, the effective Lagrangian is given by [L3] 


"E 1 
Luger = h,iv-D hy + —— | Orin + Cmag(#) Omeg | ; (3) 
2m 
where the gauge-covariant derivative D, = 0, —igsT A‘, generates the residual momentum 
ku, and 


Oxin = hy GDI hin Ong = Bh, 0.5 Gh, . (4) 


Oxin describes the kinetic energy of the heavy quark in the hadron, and Omag the heavy quark 
chromomagnetic energy. The equation of motion, iv - Dh, = 0, is exactly satisfied. The 
higher dimension operators are treated as perturbative power corrections. In the leading 
logarithmic approximation, the renormalization factor Cmag(u) is given by 


3 
—2n 
Cm (py, — roo - | 22. 7, (5) 
a(mo) 
where n; is the number of quarks lighter than the heavy quark Q. 

In the framework of HQET, it is convenient to work in the matrix representation for the 
description of the hadrons, in which wave functions of heavy hadron are only dependent on 
the heavy quark symmetry and their Lorentz transformation properties. The ground-state 
pseudo-scalar and vector heavy mesons are described by the so-called spin wave function 


pow ferJ^-03 
Mv) = Pf 7 for JS (6) 


with c" being the polarization vector. The form factors are considered as functions of the 


kinematic variable ‘ 2 2 
mg Ttm.-—q 
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Using the mass-independent normalization of the heavy meson state | H(v)) = mg *|H(P)), 
the form factors can be re-defined as 


(7(p)| uy, Q |H(v)) = 2 [fi(v - p) vu + fa(v : p) Pal; (8) 
"m 
where the dimensionless variable is p" — P .. The relation between the form factors in 
UV . 


Eqs. (M) and (B) is given by P] 


f.) = ying {Bes BP 9) 


Up MH 
2v _ 2 mj U:p 22 
fe) = Um ae ae {bite 9) tn) - <2 [i(v-p) +P fav p)] (10) 
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The form factors can be systematically analyzed order by order in terms of powers of 
1/mg. Let us begin with the expansion of the heavy-light vector current, 


(r) QI v = Fou po) 
+ Y (xp); (v) OCT). (11) 


In the limit of massless light quarks, a convenient basis of the above operators is [I4 


Ji = qy'h,, Jo = qu'h, , 
Oi-djy"iDh,, | O,— q(—iv-D)4^h,, (12) 
O;-qw'iDh,, Os = g (—iv: D) vh, 
Os = qiD"h, , Os = qd (—i D") hy. 
The corresponding Wilson coefficients are given by [4, 
1 
Bi(u)—Ci(u)—r', | Ba(w) = z Balu) = Cx(u) = 
34 4 10 16 
B —— AA d B a4 l 13 
lH) = oe — 35 — tur dn n 
28 88 20 10 
B 2 S 2 ind qoo m B 2 —9 2 —1 =. 


At the leading order of the 1/mg expansion, the matrix element of the relevant current 
qI h, can be written as 


(m(p)| aT hy |H(v)) = —Tr{ ys| Lalu - p, u) +P Lolo p, u) | T MO) }, (14) 


where the universal functions La(v - p, u) (a = a, b) depend on the kinematic variable v - p, 
but not on the heavy quark mass mg. 

At the next-to-leading order in the heavy-quark expansion, the 1/mq corrections coming 
from both the effective current and the effective Lagrangian of HQET will appear as follows. 
Matrix elements of the operators O1, O2, and O; in the effective current can be expressed 
by the generic structure 


(n(p)| q (T àD),h, |H(v)) 
= —Tr{ A Vy + Fo Pp + F3 Yp) Y + Ua Oy + F5 Pp + Fo Yp) Y f| I M(v) » (15) 


where the universal functions F;(v-p, u) (i = 1,---,6) are also mg-independent. Matrix 
elements of the operators O4, O5, and Og are not independent, and can be obtained from the 
above structures Eqs.(14) and (15). Additionally, corrections coming from insertions of the 
operators Okin and Omag into matrix elements of the leading-order currents can be described 
by six additional universal functions K,(v - p, u) (o = a, b) and S;i(v - p, u) (à = 1,---,4), 
which are defined by matrix elements of the time-ordered products 


(p) i [ay (aT (0), Oxin(y)} (9) = -Tr{ 95 (Ka + BKo)TM(v)}, (16) 


(pui | dy T (aT hy (0), Onae(y)} 18) 


= -Tef | (iS, Dal +S Oop) 5 + (iS3 Payg + Sa Cup) 5 p] IE. iid M (v) \ , (17) 
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The form factors in Eq. (B) can then be expanded as follows after doing the appropriate 
traces f], 


1 p » 
hA C La + z— L6 [Fl + (Ka Cass Se) - By(A—v-p)Lq — Bs(A—v-p) 
Q 
- 1 
x (La + Ly) + Bs| 2 — A(La + Lo) | \ + O(—), (18) 
mno 


1 - 
fo = C Li z [0 [Fh + (Ko + Cua | = By(A—v+p) Ls 


mo 
1 
+ Bo [FP eos 19] ] + OC). (19) 
MQ 
where Fi (i = 1, 3) and Sa are defined by 
Fi = F +2F; — p° Fy, F? = F — F + Fi, 
Fi = Fz + Fa +2F; — 4Fe, FÈ = Fz + F; — Fe, (20) 


Sa ES —289, + 65» + 29" $3 ; Sp = 251 = 253 + 65, . 


It should be mentioned that consequences of the heavy quark symmetry and the equations 
of motion for the heavy and light quark fields make only two (F; and Fe) of the six form 
factors F;'s independent [B]. This means that one can re-express Ff by La, Ly and Fs as 
follows 


Fi = —(A — 2w - p)La t v: pP Lr +4Fe, F? =v - pLa + ALe + 2Fe, 


- - 21 
F} = —v - p La — ALy—AFs, FR = —u-pL,— Aly —2Fe. 


Note that between the two independent universal functions from the effective current cor- 
rections only F¢(v-p, p) is relevant to the H — r matrix element in Eq. (fl) 


Fg = = (AFA): (22) 


Ill. THE LIGHT-CONE HQET SUM RULES 
Our aim is to calculate the independent form factors given in the last section by LC-HQET 
sum rules. In subsection A, the leading universal functions La (o = a, b) are calculated. 


The relevant sub-leading universal functions, Fg and Ky + CmagSa, will be calculated in 
subsections B and C, respectively. 


A. Leading order 


Let us begin with the following 2-point vacuum-pion correlation function as depicted in 
Fig. [lj 


F,(A,w) = i fatx et? (m(p)| T (aGr)nu hs), ho(0)éysd(0)} 10) 
= F (à, w) Up + F(A, w) Du » (23) 
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FIG. 1: The diagrammatic representation of the 2-point vacuum-pion correlation function Eq. (23). 
where the heavy solid line represents the heavy quark. 


where A = 2v- p and w = 2u-k. On the one hand, by inserting a complete set of intermediate 
states with the same quantum numbers as the meson between the current in the above 
vacuum-pion correlation function, we obtain the hadronic representation 


L4(À 
AI @) = 2F 711 A) ) + resonances, (24) 


where 2 = 2w - k with k being k +p. A is the heavy meson mass as defined in the HQET. 


Note that it is © that is the relevant energy of the quark system. The decay constant F is 
defined as 


(H(v)|hyPd|0) = t Furm], (25) 


with I' being an arbitrary Dirac matrix. 

On the other hand the correlation function can be calculated by expanding the T-product 
of the currents near the light-cone at z? = 0. In our calculations, the chiral limit p? = m2 = 0 
is taken. By adopting the fixed-point gauge in which 


y, A" (zr) — 0, (26) 
the following complete heavy quark propagator equals the free one, 


s) t (27) 


(fr Ch) hu(0) }10) = f datle 
The matrix elements of the nonlocal vacuum-to-pion operators can be parameterized by pion 
wave functions on the light-cone. The LC wave functions are classified according to their 
twist and the number of partons they describe. Here, let WF (n) denote the n-particle light- 
cone wave function. One can easily find that at the leading order of HQET only the WF(2)s 
are relevant. Up to twist four, the matrix elements are parameterized as follows [R, [I6]: 


(n(p)|u(z)y,ysd(0)|0) = —ifzpu [ dug P (oxlu) + z’gi(u)) 


q^puA Fi 
E = du e??? 
+f. (n, 7 i u e" P go(u) 


(tpYtryinsd(0)]0) = frin [due p(w), 


(n (p)|u(z)a,,ysd(0)|0) = j frka (Pao = pay) [ du ei Po (u) j (28) 


with ur = m2 /(m,4-m4). Here, p, is the leading twist-2 wave function, gı and g» are twist-4 
wave functions and Yp, Yo twist-3 wave functions. In the fixed-point gauge, the operator 


IIg = P exp {igs fg da z, A^ (ox)) is unity. Considering the identity 


"Yu fv = —10 py F Jw» (29) 


one can obtain FI9PT(X, ©) in terms of the light-cone wave functions: 


FLAG) = E [au [^ iate? [ween - (Fees) -o)l (30) 
Fa) = E [du [iate Zonlu) + (polu) - 20) E 
= O p^ , (31) 


where here and below u = 1 — u, and the superscript LO denotes the leading order in the 
1/mg expansion. 

Setting Fed equal F¥°®T defines the LC-HQET sum rules. The quark-hadron duality 
assumption is used to substitute the unknown resonances by the HQET result after a 
dispersion integration above some given threshold we, 


1 Im F eer A. 
resonances — — ipe aad + subtraction . (32) 
T Jwe V — wW 


The HQET spectral density can be obtained by the following double Borel transformation, 
1 PE E E 
=ImF#®7(),v) = BC? BO) Fræ Gy (33) 
T 


with the Borel transformation being defined as 


=X fixed T(n) dX” 


(34) 


whose property, B? e^*? = 5(p — 1/7), is very useful in the calculation here. Evidently, 

before doing Borel transformations, one should first perform the Wick rotation on t. 
Finally, by performing a Borel transformation Bp on both sides of the sum rule, so 

as to enhance the ground state contribution, suppress higher twist terms and remove the 


subtraction, one can obtain the sum rules, 


= We 1 
L4(X) = p | dv = ImF?"(A, v) p (35) 


GT), 


T 0 


FIG. 2: The diagrammatic representation of the 2-point vacuum-pion correlation function Eq. (4d), 
where the heavy solid lines represent the heavy quark, and the helical lines a gluon. 


where T' is the Borel parameter, and 


1 x 2 o lr 4 
ZIME (A v) = =F O(uo) | n pplu) + Ss (t) — 7392(u0) ] (36) 


1 fr 4 
Imk; (A v) = =F O (uo) U E oh (ug) = Sgi Cuo) + 9h(u0)]. — (G7) 
Here and in the following, ug = 1 — ¥, and the prime denotes derivatives. After integration 


by parts, we obtain the final sum ssl 


Lalu s p) = ze [ diu [uses (8) u HU: Pos (ü) + Zo)] e 2wvp/T 


2F 3T T 
T0 (f L T oec 
[e0 - Las] e tnr (38) 
U:p 
LEV Ur 4 2 T cesa 
L — IT Pr NN u v:p/T 
op) = gu eM ep [ du [ola o — zala) — — Gne 


+E) - 20) + LL (dn - gY(0)) | e? vw | , (39) 


where 0 = Min(1, mo) and 6 = 1 — 0. We see that the LC sum rules become meaningless 
for processes with a very soft pion which would enhance higher twist contributions. 


B. 1/mg corrections from the effective current 


To obtain the 1/mg corrections from the effective current, we consider the 2-point 
vacuum-pion correlation function containing a covariant derivative, shown in Fig. D] 


E™(,@) = if dtz e? (n(p)| T [ü(z) (T^i D), he (x), hy(0)iysd(0)} 10) 
BT (A, LM s P (A, O) Dy , (40) 
where m — 1, 3 is an index introduced for convenience, and 


(TiD), = *,yaiD? , 92D = 2D. (41) 


The hadronic representation of this 2-point correlation function can be expressed as 


T A 
qe) = 2F eU + resonances , (42) 


where the F®(A) have been defined in Eq. (RO). 

The above correlation function can be calculated directly in HQET. We choose the rele- 
vant momentum k to be parallel to v, i.e., ku = (0/2) v,, and work in the fixed-point gauge 
where A,(z) can be expressed as 


1 
Ag(z) = ae dw w Gag(wz). (43) 
0 
Eq. (#0) then can be given in terms of the above WF(2)s, and the following 3-particle wave 
functions $}, [I3]: 


(n (p) [u(x)o w 59sGae(wa)d(0)|0) 
= i fan (PaPugav — PsPugov) — (PaP- Gen — PaPrGan)| J Dorpalag or, 


(m (p) [uz )nusgsGog(w)d(0)]0) 


m Ins (gon — a] = Pa (Ipu = 2] | PoseL (oe tmo 


P (ay +wa m 
fere (pora — pata) | Dos yy (axe pat ; 


(n (p)|ü(x)o,gsGag(wx)d(0)|0) 
= de [ps (gon E a = Po( Ipu = ze) np Q1 (a,j) eter tvs) p 


px 
D ~ tay twa £ 
sr ife s Past = Pota) np Py (o;)e HIC , (44) 


where Ce = leggo G^", and Da; = da,dagda36(1 — a, — a2 — a3). The wave function 
(3x (0i) = Psr (01, 02, 03) is twist 3, and %1, Yj, #1 and ğı are twist 4. Using the identities, 


Eq. (29) and 
Yu YaB = ^IugaB — ^YouB + "189uo: — i Epai Y "rs ; (45) 


one obtains 


` z fl oo j2 T À " 
EL(A,G) = fa : du f dte SEXE ORTOS 


2 dt \ 2 
"II A E Axe 2 
+ it| Fol) — gelu) |(Bu— vu) ) — f Z Bu [es Q0) + Pau) 
sey) + it[ Pee) — go] ) ix [Ets — itai) |B 
; —iūà T 5 ES T 
—2iga(u)v, j e gt — © [pa.[ du w f dt e'2* 5t [> peel) 
- p(o) [By + 201 (a)y \ ga (46) 
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EA(,) = 4 zh du [ ated af Ary «Q0 + gilu) ] + ns) Jon 


-ia Z Ay (F[ elu) P0] + upola) ) v, 


x (v, — By) Gy (o5) e 7 7vont (47) 


where terms proportional to f? are omitted in the chiral limit. 
After performing the Ax dirus Borel transformations, we obtain the sum rules for 
the effective current corrections: 


e2A/T üe 1 
F(X) = dv = ImEZ'(A, v) e-"/T , (48) 
2F T 
where 
1 i 2 4 
-mEIQ,v) = ZO(w)( a MGE ue vee) xil (uo) ) — anuo) 
1— a 
-5 f" do, | LINE 1-a — 03,03) |, (49) 
"uo— NK. 
lmEL ) = Pn CET T ME (ug) 
" b A 0 o| 5 9x Uo; T HrPp Uo XA 0 3 Po Uo 


2 uo l-a, 
+792(uo 0) — xot uo) +f doi f dos 


0—01 


x (eyes; 1—0;-—03,03) — NE 1 — a; — 03,03) ) | , (50) 


fr 


Lm, v) = i (Geto + [nPp(uo) — at) = xen (uo) 


-i doi LL. i dos -z gen. 1 — a1 — Q3, as) |, (51) 


1 fr _ A 2 2 
—ImE(A, y) = DN | E K 3 er (uo) + MEE. E 391 (uo) ) zs xa (uo) 
l-ai 
<9 ( uo) ub sf day A 7 dos -z gen. 1 — a1 — Q3, as) |. (52) 


Note that one consequence of the equations of motion in Eq. (B1), FẸ = —F}, is explicitly 
satisfied by the sum rules. Combining the above sum rules with Eq. (22) yields, 


E 0 1— —01 
Flv- p) = -F My. | du | Pes ü )-& doy f dos 


X Qaí (01,1 — Q1 — Q3, " —— (53) 
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FIG. 3: The diagrammatic representation of the leading order vacuum-pion correlation functions 
with Okin and Omag insertions, (see Eq. BD), where the heavy solid lines represent the heavy quark, 
and the helical lines a gluon. 


C. 1/mg corrections from the effective Lagrangian 


Finally, we consider the 1/mg corrections to the correlation function Eq. (23) from the 
sub-leading operators Oy; and Omag in Eq. (H). They are depicted in Fig. B} One has 


Fy"? 0,2) = i fatzo? GET (aout), if Sy Lilu) ito) ) 0) 


GE G) 
= Fy"? (A,0)u,+ F, (Go), (54) 


where £41 = Okin + CmagOmag- By denoting 6L,/2mg as the 1/mg corrections to La, namely 
óL, = Ka + Crag Sa; (55) 
the hadronic representation of this correlator can be written as 


(u)Ha, OF OF La(A) | 2FóL,(X) 4F6AL,(A) 
da D) eg——— LIC e gas c 
(A, ©) AG 2A (2 AE + resonances, (56) 


where 0F/2mg and 0A/2mq are the 1/mg corrections for F and A [I8 [9], respectively. 
We now calculate the sub-leading correlation function Eq. (64) in HQET, starting from 
the relation 


0^ As = WF(4)s + Olas) WF(2)s, (57) 


where the symbol = indicates that this equation holds at the level of matrix elements in 
terms of LC wave functions. This relation can be found by considering Eq.(43) and the 
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equation of motion for the gluon fields, 0^G2, = =f% Ag + gshyvgT*hy. One can 
drop all the terms containing 0°Ag or A®Ag, because contributions from wave functions 
with more than 3 particles and from order a, contributions are physically quite small [B]. 
Therefore, we have 


Fy"? (Ao) = if d'z etei faty (x(p)| 7| U(x) Yyhy(x), ho(y)| - 9? + 209, A50? 


+ Cag roa | hu (y), hru(O)irsd(0) } 10) (58) 


a8 g 


Here the higher-order terms in the heavy quark propagator may be included, as displayed 
in Fig. B] (c). However, direct evaluation shows their contribution to be zero. Using the 
identities Eqs. (B9) and (49) together with 


"YuOaB = Ü (Jua YB m OuB'Ya) F Epapi YV ^fs E (59) 
"uC ag^ yv = i [Gurt gu = IJubJov) E (avus = Casus) = € uo Bv Y5 + 1 E papig” C ps ’ (60) 


and after some algebraic manipulations, we obtain 


E90, 5) — = du f ae S [S res) Po ) [B + usto Qu) 
tit | FT esu) glu) | (B vu) | ex | 2ie[ = Age (wu) — gal ) | (Gu — v) 


t 
[es (u) + gs (u) [By + usto (wo 


+ it? [2 ryo(u) — galu) (Pu — vu) ) | + (tmi) + 2igalw) 


— iiA [xai u) — itga(u) | ) By le -iuót - fpa. dss | aww [^ dte? 


À A ! À —i(l—aı —swaz)è 
x Es { 2| yy (ai) + pı (os) | — i(1 — ay — was) 3 eios) je 3)3t 


one ^ fp f dw [^ ae au — [261 (ai) — y(ai) 291 (03) | B 


+29. (c4) RD (61) 


+ 2At?gi (u) By, — ind( 


Having all necessary results at hand, we obtain the sum rules resulting from order 1/mg 
power corrections to the effective Lagrangian: 


_ 1 m, 1 Gag? (2À—v)/T x 
iL zr |, dv= Im Fa”? (A,v) + [26A/T—6F]La(A), (62) 


where the sub-leading HQET spectral density functions are given by 


(aa) T E / T /— f n = 
Emr (Av) = £e(us) [id [enit uc) + Pru) — Saf (ua) ] — 2a0[ Ha ypu) 
T / 2 / 
+ Po (uo) — xeu] }, (63) 
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Und i E D s 2 
IF," Quo) = etu) abl 362 (0) - eet) + oluo) = Sau) 


NT 7 2 6 8 12 
— 2ito| ies (uo) — Fg, (uo) + $92(uo) — xat (uo) ] + 5920) — cat (uo) 


2 6 
8 uo l-a, 1 a3 
ü a doy Ja oj | om -o 1) pilar, 1 — a1 — as, 03) 
a 
"i E) [v1 + p1 Jlor o 01 — 03,08) | | (64) 


It is interesting to note that the operator mag in fact does not contribute to the HQET 
spectral density functions even at order 1/mg. 


IV. NUMERICAL ANALYSIS 


Let us now analyze the sum rules numerically. We use the pion wave functions collected 
and carefully discussed in Ref. [B. [E]. They read, 


"OUR. fı + axi) 5 Qu -ay-ia alu) Z PIu — f — 14(u — à + 1j 
pluu) = 1+ By) Blu — a)? — 1] + BG) s [5(u — a)* — 30(u — a)? + 3], 
"no Gul m Cy) 5 Qu -ay -1]+ (Zu sg E E us} , 


5 1 6 
gılu, u) = 59 (MPU? + £009 (u) [üu(2 + 13üu) + 100? In u(2 — 3u + au) 


6 
+100” In u(2 — 3u + zu) ; 


qu p) = 5%(u)au(u =a), 
Psl u) = 360010203) 1 + ar) (Tos 3) wy 9 (i) — days — Bos F 802) 
una (u)(3a102 — 2a3 + 305) | , 
pilaf) = 309 (4) (05 — os)adl + 2e(u)(1 — 205), 
pilai H) = 1200?(u)e(u)(o1 — a2)o1a205 , 
Palisy) = 30 (iod — x) [z + 2: (4) (1 — 203]. 
Py (au, u) = ~1208%ar0205[5 + e(u)(1 — 3o3)], (65) 


where u is the renormalization scale. The scale dependence of these wave functions is given 
by perturbative QCD. 

The light-cone sum rules obtained in HQET, and the pion wave functions in Eq. (65) 
depend on the subtraction point u. However, in the lowest order of a,, the value of the 
scale u is ambiguous. Within HQET, a reasonable choice is u ~ 2A. Thus we set the 
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renormalization point at io = 1 GeV. In this case, u;(u9) = 1.65 GeV, and r(ug) = 1.08. 
The coefficients appearing in the above equation are then taken as W| 


a(o) = 0.44, — a4(uo) = 0.25, — Bo(uo) — 0.48, — Balto) = 145, 
Co(ug) = 0.10, | C4(ug) = 0.067, ó^(ug) = 0.2GeV?, elmo) = 0.5, 
fax (uo) = 0.0035 GeV?, wioluo) = —2.88, wao(uo) = 10.5, wi1(uo) ^ 0. (66) 


The pion decay constant f, = 0.132 GeV is taken from experiment. The parameters for 
heavy mesons are obtained from HQET sum rules [I3, [18, [9]: A = 0.57 GeV, F(1GeV) = 
0.46 GeV?/?, 5A = —0.35 GeV? and F = —1.92 GeV. 


A. Universal functions 


In order to make the sum rules meaningful, we must first check the existence of sum rule 
windows, which are roughly given by Agcp < T < 2A. The lower and upper limit will be 
obtained by the physical requirement that the Borel parameter 7' must be large enough to 
ensure that the higher-twist wave function contributions are suppressed, and at the same 
time small enough in order to make the resonance contribution not too large. To be specific, 
for the typical threshold values we ~ 2 GeV, we find that, setting the Borel parameter to 
0.60 GeV < T < 1.00 GeV for the sum rules La and Ly, the twist-4 wave functions give 
contributions less than 1796 and 396, while the resonance contributions are lower than 2796 
and 6%, respectively. The two sum rules are quite stable in this region, as shown in Fig. H(a) 
and (b) where v - p is fixed at 2.0 GeV. Nevertheless, our numerical calculations show that 
the stability of the sum rules is not very sensitive to a change of we. In order to be consistent 
with HQET sum rule calculations, we adopt the same stability region of we as that from the 
A sum rule [I3, [8]. i. e., we = 2.0 + 0.3 GeV. Fig. Df(a) and (b) present the final results for 
Lalu- p) and L,(v- p), respectively, where the central value of the Borel parameter T = 0.80 
GeV is used. Note that the LC-HQET sum rules are meaningless in the soft pion region as 
mentioned above. To stay away from the soft pion region, we take v - p > 1.2 GeV. 

For the convenience of further applications, we parametrize the results by the following 
formulae in the region 1.2 GeV < v: p < 2.64 GeV, 


1 
LE S = 67 
a p) EN TANN (67) 


b 
LE* (v -p) = bo thu: p+ TP (68) 


Best fit values in which the maximal errors are less than 296 yields the following values of 
the parameters: 


We [GeV] lao [GeV-?] ay [Gev-3/2] az [GeV 9/2] bo [GeV] bı [Gev-¥/2] bə [Gev3/2] 
2.3 -0.78 4.00 0.891 0.403 -0.0221 -0.0385 
2.0 -3.03 6.71 0.286 0.409 -0.0263 -0.0448 
1.7 -6.25 10.9 -0.743 0.320 -0.0049 0.0343 


Next we consider the 1/mg corrections from the effective current. In the relevant universal 
function Fy, only two twist-3 wave functions appear. Thus, in order to determine the 
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FIG. 4: Variation of the sum rule predictions on La,» with the Borel parameter T at various values 
of the continuum threshold parameter we. 


sum rule window, especially its lower limit, one may consult the sum rules for Fl, that 
will contribute to the form factors directly while the renormalization-group effects are not 
considered. For threshold values of we ~ 1.8 GeV, requiring that both the twist-4 wave 
functions and the resonance contributions do not exceed 40%, the sum rule window can 
be determined to be 0.50 GeV < T < 0.70 GeV, in which we find that the resonance 
contributions for Fg are less than 8%. The numerical results for the universal function F¢ as 
a function of the Borel parameter T are shown in Fig. f|, where v - p is fixed to be 2.0 GeV. 
Fig. [/ presents the results for Fs(v- p), where the central value T = 0.60 GeV is taken for 
the Borel parameter. 

The results can be parametrized by the formula below valid in the region 1.2 GeV < 
v- p « 2.64 GeV, with maximal errors less than 3%, 


; 1 
pl] c" = === a, 69 
6 (v p) C9 +c v- p+ c (v - p)? ( ) 
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FIG. 5: The universal functions La (v : p) obtained from LC-HQET sum rules at various values 
of the continuum threshold parameter we. 


The parameters c; are given by 


We [GeV] Co [GeV-3/?] Cy [GeV - 9/7] Co [GeV-7/?] 
2.1 -112 30.5 -4.93 
1.8 -86.1 1.54 2.22 
1.5 -57.6 -34.0 11.7 


We next analyze the sum rules resulting from the 1/mg insertions due to the power 
corrections to the effective Lagrangian. We find that, for the sum rule for à L,, the twist-4 
wave function contributions are very small everywhere in the physically appropiate region 
of T, but the resonance contributions grow rapidly as T becomes large. Requiring the latter 
contributions to be less than 40% yields T < 0.55 GeV. Consequently, we set the range of 
the Borel parameter at 0.35 GeV < T < 0.55 GeV. The results of the sum rule for dL, are 
very similar to the leading order sum rules. In the sum rule window 0.60 GeV < T < 1.00 
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FIG. 6: Dependence of the sum rule predictions of Fg on the Borel parameter T at various values 
of the continuum threshold parameter we. 
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FIG. 7: The sub-leading universal function Fg(v - p) obtained from the LC-HQET sum rules at 
various values of the continuum threshold parameter we. 


GeV and for threshold values we ~ 1.8 GeV, the twist-4 wave function contributions are less 
than 11% and the resonance contribution does not exceed 18%. The stability of the sum 
rules for dL, (a = a, b) with regard to variations of the Borel parameter T are shown in 
Fig. B(a) and (b), where v - p is fixed at 2.0 GeV. Fig. B(a) and (b) present our final results 
for ÓL,(v- p) and 6L,(v-p), respectively, where the corresponding central value of the Borel 
parameter T' — 0.45 and 0.80 GeV are used. It should be mentioned that, in these two 
sub-leading universal functions, the sub-leading HQET spectral density functions give small 
and negative contributions to the form factors, while the corrections coming from the ôF 
and 6A contributions are large, and positive when the sum is taken (for ôLa, the proporation 
is about —1 : 2, and ôL, about —1 : 5 ). 

We have parametrized these results by the following formulae in the region 1.2 GeV « 


Lum dms 
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FIG. 8: Dependence of the sum rule predictions La,» on the Borel parameter T at various values 


of the continuum threshold parameter we. 


v-p< 2.64 GeV, 


| 1 
óL (v.p) = 


"E b 
óLP*(v.p) = borho p ——. 

U-p 
The best fitting parameters are 


We [GeV] |đo [GeV-?/72] G1 [GeV-9/7] Gq [Gev- 7/2] bo [GeV3/2] 
21 -4.56 10.6 -1.59 0.319 
1.8 1.62 3.57 0.088 0.129 
1.5 6.05 -2.16 1.55 0.539 
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FIG. 9: The sub-leading universal functions 9L, (v - p) obtained from LC-HQET sum rules at 
various values of the continuum threshold parameter we. 


where the maximal fitting errors are less than 5%. 

Finally we remark on the errors in the sum rules. The uncertainties of the continuum 
threshold parameter we induce errors (from both the LC-HQET sum rules themselves and 
the sum rules of A and F [L3]) less than + 10% in the final sum rule results. However, one 
should keep in mind that the sum rule method typically has a +(10 ~ 30)% uncertainty 
resulting from the duality assumption, uncertainties in the wave functions and other input 
parameters. 


B. Form factors to order 1/mg 
Having obtained the leading and the relevant sub-leading universal functions, we can 


construct the form factors fi,»(v- p) from Eq. (18) for the B — 7 transition to order 1/m;. 
The results are shown in Fig. [LO] (a) and (b), where we take m, = 4.7 GeV. For comparison, 
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the leading order results corresponding to we = 2.0 GeV are also shown in this figure. Here 
the sub-leading universal functions roughly give a 11 — 2396 and 5 ~ 12% enhancement for 


JPT (v - p), respectively. 


the form factors fP "(v - p) and 


ph > T(y.p) Gev!2 
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FIG. 10: The form factors fE (v : p) obtained from the LC-HQET sum rules at various values 
of the continuum threshold parameter we. 


V. COMPARISON WITH OTHER MODEL CALCULATIONS, DISCUSSION 
AND SUMMARY 


It is interesting to compare our results with those of other approaches. First let us look at 


the soft pion limit. The heavy meson chiral perturbation theory (HMyPT) [B0] (see also P1]) 
describes fi(v- p) and fa(v- p) in a single pole form in the soft pion region. Moreover, the 
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1/mge corrections do not change this behaviours [Bl 22]. In the case of B — 7 transitions, 
they can be expressed by 


C,F oF s 

Bor 1 ÉL gv'P 

Li -— 1 1 Se 2 

Ree p) 27. ( fo) | 0) p+ Ag f M 
CF OF à 

Bor . = 1 _ IVP, 

2 HMxPT (v p) E 2 fr (1 2my ) v:pdcAg' e) 


where Ag ~ Mp — mg = 0.045 GeV, g ~ 0.3 is the coupling of the pion to the heavy 
meson. It is easy to see that our result for fP^7(v-p) does not match to the HMyPT result. 
Concerning f?~*(v - p) our sum rule result Eq. (89) can be seen to reasonably well match 
on to the soft pion result considering the uncertainties in the values of F and g. By taking 
g = 0.28, one can find that the extrapolation of HMyP to large v: p matches quite well with 
the sum rule calculations at intermediate pion energies as shown in Fig [L| (a). For practical 
purposes, we use the following Gaussian-type function to make a smooth connection between 
the sum rule result of fP^"(v. p) and that of HMyPT, 


f») = S (1+ 2) [go ge tmm]. (74) 


(for 0.25 GeV < v-p< 1.2 GeV) 


where go = 0.125, g; = 0.751, go = 3.0 GeV? and g = 0.23 GeV were obtained by 
matching both sides. This is plotted in Fig. (b). It should be mentioned that the 
matching results obtained here correspond to the case of central values of we. For further 
support of our matching procedure we have included in the lattice NRQCD results from the 
JLQCD collaboration [23]. They can be seen to be well consistent with our matching result. 

The large pion energy limit is another interesting limit to consider because the large 
energy effective theory (LEET) [24 B3] provides some model-independent information on 
this kinematical region. In the B — 7 semileptonic decay, the light non-spectator quark gets 
a large amount of energy from the decay of the bottom quark. The light spectator quark 
system interacts with the energetic quark mainly at the energy scale Aocp which is basically 
fixed by the size of hadrons. New symmetries appear in the limit of v - p/Aocp — oo. They 
are subject to corrections by hard gluon exchange [6]. In the leading order of the heavy 
quark expansion our results are 


E 


La(oo) = 0, Ly(oo) E ye) lo, (75) 
where i w 
I; = ze | drr'e™” , (76) 
4 0 
This agrees with LEET. For the 1/mg corrections, our results are 
fr Hr 
Fo) = - Ee) ho, 
óLs(oo) = 0, (77) 


óLy(oo) = i T | (26A/T — ôF) I -T 13) | &«(1). 
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FIG. 11: The form factors JE 3" (v-p) as obtained from LC-HQET sum rules including a matching 
to the HMxPT result in the soft pion region. 


This can be compared with the LEET results only after the hard gluon effects have been 
incorporated into these and can thus only been done in the future. There even is a school 
which assumes that the perturbative contribution is dominant for the B — 7 transition 28]. 
It is obvious that the LEET limit deserves further studies which certainly will be done in 
the near future. 

It is also interesting to compare our method and our results to the full QCD LC sum 
rule calculation. The authors of Ref. [29] in fact considered the heavy quark limit of their 
full QCD LC sum rule calculation. Their results differ from L,(v- p) and Li(v- p) only by 
a simple transformation which were given in Eqs. (47) and (46) in their paper. There is a 
subtle difference from our results, though, which can be seen by letting x(1 — u)/u = v. It 
is only after taking v/v-p — 0 that the results of the two calculations fully agree. This limit 
may need further understanding. Numerically the QCD LC sum rules gave stable results for 
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q? < 17 GeV? which is consistent with the HQET LC sum rules requiring v - p > «4/2 ~ 1 
GeV. These two methods could be in principle the same, provided that all the sub-leading 
corrections have been included in the calculations. Lacking such powerful calculations, the 
effective theory calculates a physical quantity in a most thorough and least complicated 
way through clearly separating the perturbative and nonperturbative parts of the quantity. 
Therefore in certain appropriate regions of the pion's phase space the HQET calculation 
may give more reliable results. Practically speaking the method is to calculate form factors 
at the hadronic scale, and the perturbative contribution is accounted for by multiplying in 
renormalization factors. 

It is physically useful to reconstruct the conventional form factors defined in Eq. (|l] 
from the results combining LC-HQET sum rules and HMxPT, by using the relation Eq.(B). 
In Fig. (a) we present f?~*(q?), which is directly measurable in semi-leptonic decay 
involving light leptons. In Fig. (b) we present our result for the scalar form factor 

(77(g?). The scalar form factor contributes to the decay B — «7v, and also enters the 
factorized amplitudes in non-leptonic two-body B decays. For comparison, the results from 
the lattice QCD simulations by APE BO, UKQCD Bi], FNAL and JLQCD are 
also shown in the figures. 

In Table [| we compare our results to the results of other model calculations which includes 
the conventional HQET QCD sum rule calculation of Ref. [B3]. We find that our results 
for fP^"(g?) are about 14 ~ 33% larger than those from the full QCD LC sum rule [I] 
B4] calculation which also includes a, corrections. 


TABLE I: Comparison between different approaches for the universal functions La(v-p) and Lp(v-p), 


and the form factor fP—7(q?) at some different kinematical points. 


L«(v:p)[GeV!2] Lelu: p) [GeV/7] f«(q?) 
v : p [GeV] or q? [GeV?] 1.50 2.64 1.50 2.64 0.0 6.0 12.0 
LC-HQET SR (NLO) — — 0.42 +0.04 054 0.77 
LC-HQET SR (LO) 0.13 0.06 0.34 0.32 0.36 0.47 0.68 
HQET SR in LO 0.15 0.13 0.25 0.21 0.24 0.32 — 
LC-QCD SR — — 0.28 +0.05 0.42 0.66 


Finally, in order to be concrete, we give our prediction for the decay width of the decay 
B? — r-e*vy, using our form factors. We obtain 
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Vu (1.46 + 0.30) x 1079 GeV. (78) 


~ [4.08 x 10-3 


By taking |Vaæ| = (4.08 + 1.18) x 10? from inclusive measurement of B — X,tv [B5], 
we get T = (1.46 + 0.307022) x 10776 GeV, where the uncertainties are from that of form 
factors and |V,;| respectively. On the other hand, from the experimental result given in 
Ref. BG], Tso = (1.55 + 0.03) ps^! and Br(B? — «^ e*v,) = (1.8 E 0.6) x 1074, we extract 
Viv] = (2.947033 + 0.50) x 1073 with the uncertainties from the form factors and the 
experiments, respectively. 

Let us also comment on Ref. which gave leading order results for La. In addition 
to that the analytical expressions of the sum rules differ from ours by a factor of 2, the 
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FIG. 12: The form factors Jea Ca from combining LC-HQET sum rules and HMyPT, and from 
the lattice QCD simulations by APE [Bd], UKQCD [Bi], FNAL [BJ] and JLQCD 23}. 


numerical input used in Ref. [[1] is quite different, such as the choice of the energy scale u 
and the determination of the sum rule window. This naturally affects the final results on 
Lalv- p). 

To summarize, we have applied the LC sum rule method to calculate the B — mlv 
weak decay form factors to order 1/mg in the framework of HQET. We have calculated the 
leading and the relevant sub-leading universal form factors. Our form factor results have 
been matched to the appropiate soft pion results. We have also discussed the large pion 
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energy limit of our results. The full QCD LC sum rules and lattice QCD results have been 
compared with. In the future we are planning to include perturbative QCD corrections, 
try to incorporate the large energy effective theory into the QCD sum rule technique and 
perform detailed phenomenological analysis. 


Acknowledgments 


We would like to thank Chao-Shang Huang, Stefan Groote and M. Beneke for helpful 
discussions. C.L. acknowledges support from the Alexander von Humboldt Foundation. 
This work was supported in part by the National Natural Science Foundation of China with 
the grant no. 10075068, and by the BEPC National Lab Opening Project. We would like 
to thank ICTP for hospitality where part of this work was done. 


1] For reviews, see V.M. Braun, hep-ph/9801222; A. Khodjamirian and R. Rückl, 
ph/9801433. 
For a review, see M.B. Wise, lhep-ph/9805468. 


M. Beneke, G. Buchalla, M. Neubert and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999). 
Y.Y. Keum, H.-n. Li and A.I. Sanda, Phys. Rev. D 63, 054008 (2001). 

G. Kramer, G.A. Schuler and T. Mannel, Z. Phys. C 51, 649 (1991). 

G. Burdman, Z. Ligeti, M. Neubert and Y. Nir, Phys. Rev. D 49, 2331 (1994). 

C.-S. Huang, C. Liu, and C.-T. Yan, Phys. Rev. D 62, 054019 (2000). 

V.M. Belyaev, V.M. Braun, A. Khodjamirian and R. Rückl, Phys. Rev. D 51, 6177 (1995). 
V.M. Belyaev, A. Khodjamirian and R. Rückl, Z. Phys. C 60, 349 (1993); 

A. Khodjamirian, R. Rückl, S. Weinzierl and O. Yakovlev, Phys. Lett. B410, 275 (1997); 
E. Bagan, P. Ball and V.M. Braun, Phys. Lett. B 417, 154 (1998); 

T. Huang, Z. Li and X. Wu, Phys. Rev. D 63, 094001 (2001). 

10] Y.-B. Dai and S.-L. Zhu, Phys. Rev. D 58, 174009 (1998). 

11] W.-Y. Wang and Y.-L. Wu, Phys. Lett. B 515, 57 (2001). 

12] A. Ali, V.M. Braun and H. Simma, Z. Phys. C 63, 437 (1994). 

13] For a review, see M. Neubert, Phys. Rep. 245, 259 (1994). 

14] A.F. Falk, M. Neubert and M.E. Luke, Nucl. Phys. B 388, 363 (1991). 

15] A.F. Falk and B. Grinstein, Phys. Lett. B 247, 406 (1990). 

16] V.M. Braun and LB. Filyanov, Z. Phys. C 48 239 (1990); 

P. Ball, JHEP 9901, 010 (1999). 

17] Y.-B. Dai, C.-S. Huang, C. Liu and C.D. Lü, Phys. Lett. B 371, 99 (1996); 

Y.-B. Dai, C.-S. Huang, M.-Q. Huang and C. Liu, Phys. Lett. B 387, 379 (1996). 

18] M. Neubert, Phys. Rev. D 45, 2451 (1992); 

E. Bagan, P.B. Ball, V.M. Braun and H.G. Dosch, Phys. Lett. B 278, 457 (1992). 

19] M. Neubert, Phys. Rev. D 46, 1076 (1992); Phys. Lett. B 389, 727 (1996); 

P. Ball, Nucl. Phys. B 421, (1994) 593; 

20] M.B. Wise, Phys. Rev. D 45, R2188 (1992); 

T.-M. Yan et al., Phys. Rev. D 46, 1148 (1992), Erratum: Phys. Rev. D 55, 5851 (1997); 
G. Burdman and J.F. Donoghue, Phys. Lett. B 280, 287 (1992). 


O OND oO FB Ww bd 


25 


21 
22 
23 
24 
25 


26 


27 
28 


29 
30 
31 
32 
33 
34 


35 


36 


C.A. Dominguez, J.G. Korner and K. Schilcher, Phys. Lett. B 248, 399 (1990) 

C.G. Boyd and B. Grinstein, Nucl. Phys. B 442, 205 (1995). 

S. Aoki, et al, Phys. Rev. D 64, 114505 (2001). 

M.J. Dugan and B. Grinstein, Phys. Lett. B255, 583 (1991). 

J. Charles, A. Le Yaouanc, L. Oliver, O. Péne and J.-C. Raynal, Phys. Rev. D 60, 014001 
(1998). 

M. Beneke and Th. Feldmann, Nucl. Phys. B 592, 3 (2001); 

C.W. Bauer, S. Fleming, D. Pirjol and I.W. Stewart, Phys. Rev. D 63, 114020 (2001). 

B. Grinstein and D. Pirjol, Phys. Lett. B 533, 8 (2002). 

R. Akhoury, G. Sterman and Y.-P. Yao, Phys. Rev. D 50, 358 (1994); H.-n. Li and H.-L. Yu, 
Phys. Rev. Lett. 74, 4388 (1995); T. Kurimoto, H.-n. Li and A.I. Sanda, Phys. Rev. D 65, 
014007 (2002); Z.-T. Wei and M.-Z. Yang, 

A. Khodjamirian, R. Rückl and C.W. Winhart, Phys. Rev. D 58, 054013 (1998). 

A. Abada et al, Nucl. Pyhs. B 619, 565 (2001). 

UKQCD Collaboration, Phys. Lett. B 486, 111 (2000). 

A. El-Khadra et al, Phys. Rev. D 64, 014502 (2001). 

P. Colangelo and P. Santorelli, Phys. Lett. B 327, 123 (1994). 

A. Khodjamirian, R. Rückl, S. Weinzierl, C. W. Winhart and O. Yakovlev, Phys. Rev. D 62, 
114002 (2000). 

A. Bornheim, et al. (CLEO Collaboration), Phys. Rev. Lett. 88, 231803 (2002); 

For a recent review, see H. Schwarthoff, lhep-ex/0205015 

J.P. Alexander, et al. (CLEO Collaboration), Phys. Rev. Lett. 77, 5000 (1996); 

D.E. Groom, et al. (Particle Data Group), Eur. Phys. J. C 15, 1 (2000). 


26 


